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We analyze the proper inclusion of electric-field-induced spin splittings in the framework of the 
envelope function approximation. We argue that the Rashba effect should be included in the form 
of a macroscopic potential as diagonal terms in a multiband approach rather than the commonly 
used Rashba term dependent on k and electric field. It is pointed out that the expectation value of 
the electric field in a subband is sometimes not unique because the expectation values can even have 
opposite signs for the spin-split subband components. Symmetric quantum wells with Dresselhaus 
terms and the influence of the interfaces on the spin splitting are also discussed. We apply a well 
established multiband approach to wide modulation-doped InGaSb quantum wells with strong built- 
in electric fields in the interface regions. We demonstrate an efficient mechanism for switching on 
and off the Rashba splitting with an electric field being an order of magnitude smaller than the local 
built-in field that determines the Rashba splitting. The implications of our findings for spintronic 
devices, in particular the Datta-Das spin transistor and proposed modifications of it, are discussed. 
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I. INTRODUCTION 

The interaction of a spin with a magnetic field or mag- 
netic ions has been studied extensively over the years. A 
surprisingly efficient mechanism to produce a spin split- 
ting is utilizing the combination of inversion asymmetry 
and spin-orbit interaction, A One convenient way to con- 
trol the spin splitting is the Rashba effect which results 
from structure inversion asymmetry (SIA)j2 An applied 
or built-in macroscopic electric field is seen in the frame of 
a moving electron as having a magnetic- field component 
and yields a spin splitting. In this way one can utilize 
many of the mechanisms from conventional electronics 
which is mainly controlled by electric fields. The idea of 
using the spin of the carriers in addition to its charge has 
resulted in a research area called spintronicsi^ Another 
spin splitting mechanism resulting from the lack of inver- 
sion symmetry of the zinc-blende lattice [bulk inversion 
asymmetry (BIA)] is called the Dresselhaus effect. 4 Both 
the Rashba and the Dresselhaus effects are frequently in- 
cluded via terms linear in the wave vector k. However, 
they are the lowest-order terms of more accurate expres- 
sions that are obtained from multiband envelope function 
theory* 1 - 

In Sec. II we will recapitulate the foundations of the 
commonly used envelope function approximation^ in or- 
der to set the ground for an analysis of the proper in- 
clusion of SIA and BIA within the framework of this 
approximation. We make here the important distinc- 
tion between slowly and rapidly varying potentials. In 
Sec. Ill we apply the multiband theory to an interest- 
ing system, wide n-type modulation-doped (MD) quan- 
tum wells (QWs) with strongly nonuniform electric fields. 
Here there are strong built-in electric fields of opposite 
signs at the two interfaces and one can expect a strong 
Rashba effect. It has frequently been assumed that one 
needs a strong applied bias to get a substantial Rashba 



splitting but interesting things can happen also for small 
or moderate bias. We show how the built-in electric 
fields in modulation-doped quantum wells can be uti- 
lized while applying a much smaller external field. For 
very small asymmetry interesting anticrossing phenom- 
ena occur. Furthermore, the spin splitting due to the 
Dresselhaus effect in symmetric quantum wells is found 
to be qualitatively different in modulation-doped quan- 
tum wells compared to square wells. In Sec. IV we dis- 
cuss the implications of our results for spintronic devices, 
in particular the Datta-Das spin field effect transistor. 7 
It has for a long time been considered as a prototype of 
a spintronic device but unfortunately the efforts to im- 
plement it in practice have not yet been very successful. 
Finally, in Sec. V we discuss the results and conclude. 



II. THEORY 

The envelope function approximation (EFA) has been 
widely used during several decades. Under the name 
effective-mass theory, it was first applied to shallow- 
impurity states in bulk semiconductors 4^ The starting 
point is that the problem with the band structure in the 
pure bulk material is assumed to have been solved. Ac- 
cording to Bloch's theorem the total wave function for 
band n is given by 

^nk(r) = e ik - r Unk (r) (1) 

where u„k(r) has the periodicity of the lattice. We then 
introduce a perturbation U (r) . An essential assumption 
in the derivation is that it should be slowly varying on the 
scale of the lattice constant. This assumption does not 
always hold in the cases where the EFA has been applied. 
The advantage of the EFA is that one can avoid the ex- 
plicit inclusion of the cell periodic potential. Only the 
slowly varying perturbation U(r) enters a Schrodinger- 
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like equation. With the perturbation the total wave func- 
tion can be expanded 

V(r) = X>( r Kk(r) (2) 

n 

where the summation is over all the energy bands. / n (r) 
is called an envelope function and the EFA gives a simple 
prescription for the effective Hamiltonian operating on 
the envelope function /„ only,— £ 

H = E{-iV) + U{v) (3) 

The kinetic-energy operator E is obtained from the 
bulk band structure E(k). In general it is a matrix whose 
eigenvalues gives the energy-band dispersion in the bulk. 
For the perturbed problem k is replaced by — iV, where 
V is the gradient operator. In the case of a quantum well 
grown along the z direction, it is sufficient to replace k z 
by the operator —id/dz while k x and k y remain good 
quantum numbers. The perturbation potential U(r) is 
added along the diagonal of the matrix. 

For the impurity case the potential is not slowly vary- 
ing in the unit cell containing the impurity. For s states 
having a finite amplitude at the origin this sometimes 
gives important deviations (often called central-cell cor- 
rections) from the predictions of the EFA in its simplest 
form»iS The ground state for a donor in Si is, for exam- 
ple, split up into three levels where the energy separation 
between these levels can be comparable to the predicted 
ground state energy. For p and d states in Si the EFA 
works excellently and it also works well for the ground 
state in direct-gap semiconductors ^ 

The EFA has also been used frequently for quantum 
well heterostructures. Here the potential changes rapidly 
near the interface between two materials. The range of 
this potential change is of the order of the lattice con- 
stant. Thus, at first sight, it appears that the EFA would 
not be applicable. In spite of this, the agreement between 
its predictions and experimental results has turned out 
to be quite good if the EFA is applied properly with- 
out unnecessary approximations. The reason for this 
was examined by Burfei^ in a series of papers and led 
to a new set of boundary condition s 13 ' 14 nowadays called 
the Burt-Foreman boundary conditions. These bound- 
ary conditions were first derived for the 6x6 Hamilto- 
nian describing the valence bands^ 3 - and later a prescrip- 
tion for extending it to the 8x8 Hamiltonian including 
the conduction band was givenJ^ Previously it was com- 
mon to symmetrize the order between operators and spa- 
tially varying parameters but the Burt-Foreman bound- 
ary conditions were derived with consideration of the cell- 
periodic wave functions w„k(r). Meney et al^ analyzed 
different envelope function approaches and found support 
for the Burt-Foreman boundary conditions. Burt's anal- 
ysis also explained why the EFA works quite well even 
for narrow quantum wells. However, it should be kept in 
mind that caution is necessary when applying the EFA to 
interface regions. In the EFA the interfaces have usually 
been taken as abrupt steps. A rapidly but continuously 
varying potential has also been considered by Stern and 
Das Sarma 1 ^ but the influence on the subband energies 
was found to be quite small. 



The summation in Eq. (2) should in principle be over 
all the bands. In practice one selects a finite number of 
important bands whose interaction is included exactly 
in the matrix while the "remote" bands are included 
perturbativelyJ^ A larger number of bands included in 
the matrix gives an accurate description in a larger k 
range. A common choice that we apply in this paper 
is to include the conduction, heavy-hole, light-hole, and 
split-off bands in an 8 x 8 matrix. For a symmetric struc- 
ture this includes a twofold spin degeneracy. 

For even more accurate results, 14 x 14 and 16 x 16 
matrices have been considered^ Wissinger et alJ^ have 
performed calculations for an asymmetric GaAs quan- 
tum well both with the 8x8 and the 14 x 14 matrices, 
and compared to Raman-scattering experiments^ The 
deviation between the two models was rather small and 
comparable to the deviation from the experimental re- 
sults. In the present case we consider InGaSb with a 
small band gap where the 8x8 matrix is expected to be 
a good approximation. On the other hand, it can be suf- 
ficient to use 6x6 matrices in which the conduction band 
or the split-off band is among the remote bands. The fre- 
quently applied Luttinger-Kohn Hamiltonian^ includes 
the heavy-hole and light-hole bands in a 4 x 4 matrix. For 
a description of electron subbands, it is convenient to use 
a one-band (two-component) approximation in which all 
the other bands are considered as remote and included 
via a modification of the free-electron mass to an effec- 
tive mass. As we will see below, inclusion of spin effects 
in a one-band model leads to some complications. 

The recapitulation above of the essence of the EFA 
serves as the basis for analyzing the inclusion of 
asymmetry-induced spin phenomena. The inclusion of 
the Dresselhaus effect 4 - seems clear. It stems from the 
microscopic structure of the bulk material, influences the 
cell-periodic part w„k of the wave function, and results 
in a modified bulk band structure. Thus it is appropriate 
to include it as ^-dependent terms in the matrix which, 
after the replacement k — > — iV, becomes the kinetic- 
energy operator. Several terms of different order in k 
enter the 8x8 Hamiltonian^ For electron subbands the 
lowest-order term is linear in k: 

H D = (3{k x (J x - kyOy), (4) 

where a x and a y are Pauli matrices, and (5 is a material 
constant giving the strength of the Dresselhaus effect. 

The inclusion of the Rashba effect^ is less uncontrover- 
sial. It stems from a slowly varying macroscopic electric 
field and, according to the principles of the EFA, it should 
be included as a z-dependent potential along the diago- 
nal of a matrix of sufficient size. The Rashba effect in 
p-channel Si metal-oxide-semiconductor field-effect tran- 
sistor (MOSFET) structures was already included in this 
way in the 1970s and good agreement with experiment 
was found^ Using a multiband matrix for the kinetic en- 
ergy, the inclusion of an asymmetric potential results in 
a spin splitting for finite values of the in-plane wave vec- 
tor without inclusion of any special ^-dependent terms 
[cf. Eq. (5) below]. The spin-orbit interaction is im- 
plicitly included via the coefficients of the fc-dependent 
elements in the matrix. They contain matrix elements of 
the spin-orbit interaction with respect to the cell-periodic 
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wave functions u„k and can be evaluated theoretically. 4 
However, a higher accuracy can often be obtained from 
cyclotron resonance experiments,— and in practice ex- 
perimentally determined effective masses and Luttinger 
parameters^ are inserted if they are available. 

Cyclotron resonance experiments for a two- 
dimensional hole gas at a modulation-doped 
GaAs/AlGaAs interface were performed by Stormer et 
al£^ The roughly triangular potential with a strong 
Rashba effect gave rise to two clearly different masses 
ascribed to the two components of the spin-split heavy- 
hole subband. These results were in very good agreement 
with calculations with the Luttingcr-Kohn Hamilto- 
nian, where the measurable energies of the allowed 
Landau-level transitions were explicitly calculated^ 4 - 

An alternative to the numerical solution of a multi- 
band problem with suitable boundary conditions is a 
process called "downfolding." Starting from a multiband 
matrix, one can derive various 2x2 Hamiltonians . 25 ' 26 ' 27 
The commonly used Rashba term can be derived as 
the lowest-order term»i For the electric field in the z- 
direction, it becomes 2 - 

H R = a{k x a v - k y a x ) (5) 

Here a is often used as an input parameter taken from 
experiment. Several experiments (see, e.g., Refs. I20I. l28l. 
and I29I ) have aimed at determining this Rashba coeffi- 
cient for different materials. Various theoretical expres- 
sions have been derived for a. A simple expression that 
we will use as a reference is 2 ^ 

h 2 eA(2E g + A) 
a 2mE g {E g + A)(3E g + 2A) ' U 

where < e > is the expectation value of electric field in 
the quantum well and the barriers. 

The Rashba term is not really consistent with the prin- 
ciples of the EFA. It is a kind of hybrid including both 
the wave vector k and the potential. The problem arises 
from the fact that the s-like conduction band gets its 
spin-orbit coupling from the interaction with the valence 
bands, which are included in this approximation among 
the "remote" bands. A further problem is that Eq. (6) 
implicitly assumes that the electron subband has a well- 
defined expectation value but as will be discussed below, 
the two components of the spin-split subband (hence- 
forth denoted spin subbands) can have clearly different 
expectation values. 

A special problem is how the interfaces of a quantum 
well should be included. In an asymmetric quantum well 
the penetration of the wave function into the left and 
right barriers becomes different. At a first glance it seems 
natural to treat the complete conduction-band profile as 
the relevant potential. Zawadzki and Pfeffer— have in- 
cluded the conduction band offsets in < e > [cf. Eq. (6)] 
and denoted it the " average electric field." Using the fact 
that no force acts on a bound state it has been arguecU^ 
that the contribution from the interfaces would largely 
cancel that of the electric field in the quantum well. How- 
ever, in this respect it is important to distinguish be- 
tween the total wave function and the envelope functions 
for which this argument does not necessarily hold. The 
different effective masses in well and barrier, and spin 



dependent boundary conditions make this " average elec- 
tric field" yield a non-zero but small contribution to the 
Rashba spin splitting. 

Lassnig22 argued that the valence-band profile, includ- 
ing the band offsets, determines the Rashba effect for 
conduction electrons. This gives an interface contribu- 
tion of the same sign as that of the electric field in the 
quantum well and barriers. This interface contribution 
has been evaluated analytically by downfolding the other 
bands on the conduction band and resulted in matrix ele- 
ments of the step-like valence-band edge3£ However, this 
approach of treating the interfaces like infinite electric 
fields seems dubious against the discussion above, where 
it was pointed out that the interfaces are a weak point in 
the EFA. A crucial nontrivial factor in the downfolding 
procedure is how the order between differential opera- 
tors and spatially varying material parameters should be 
chosen to be compatible with the Burt-Foreman bound- 
ary conditions . 13 ' 14 So far downfolding approaches have 
usually used the ad hoc operator symmetrization that 
many workers have abandoned in multiband approaches. 
We will analyze the interface contribution below in our 
multiband approach using Eqs. (5) and (6). 

This interface contribution is fundamentally different 
from another interface contribution where the micro- 
scopic structure at the interface is taken into account as 
an additional source of inversion asymmetry. The proper 
inclusion of this short-range potential in the framework 
of envelope function theory is not trivial and has been 
subject to debated The effect is particularly strong in 
"no-common-atom interfaces" such as InAs/GaSb where 
the two constituents have no atom in common^ We have 
neglected it in the present calculations. 

Equations (5) and (6) predict a linear increase in the 
Rashba spin splitting with the in-plane wave vector k. 
Numerical result s 25 ' 32 give a Rashba splitting that is non- 
linear and can even be nonmonotonic. Yang and Chang 32 
have recently published numerical and approximate an- 
alytical solutions to the Rashba effect. By inserting the 
actual subband energy e(k) instead of the bulk band- 
edge energies as in simpler models, their analytical model 
qualitatively reproduces the nonlinear k dependence. 
However, the analytical model overestimates the spin 
splitting. Their numerical results for an Ino.53Gao.47As 
quantum well between Ino.52Alo.4sAs barriers agree very 
well with what we obtain in our approach. The quantita- 
tive agreement for Hgo.74Cdo.26Te/HgTe structures is less 
good but the qualitative behavior of an increase followed 
by a decrease in the spin splitting is reproduced. We find 
a similar nonmonotonic behavior for an InGaSb/InAlSb 
structure with the same composition as in this paper but 
with a uniform electric field. On the other hand, for a 
GaAs/AlGaAs quantum well we find a monotonic but 
clearly nonlinear behavior, in agreement with Ref. [25l . 
This is consistent with the explanation 3 - 2 - that the non- 
monotonic k dependence is caused by the increased en- 
ergy separation to the light-hole and split-off bands, an 
effect that is less pronounced in GaAs with a fairly large 
band gap. It can be noted that the dependence of AE 
on k is similar to the dependence for holes in Ref. [H of 
Ak at the optimal energy on the electric field which was 
explained in a similar way. 
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We have started from the properties of the bulk mate- 
rials and, via the boundary conditions, obtained the sub- 
band dispersions. An alternative approach is to solve the 
subband problem for k — and expand the wave func- 
tions for finite k in eigenfunctions for k = (the mim- 
ic • p method).— This approach can be illuminating also 
for spin phenomena.! For a modulation-doped interface 
it was found that, however, the convergence with respect 
to number of basis functions was surprisingly slow^ 



III. RESULTS 

We here consider an interesting system for studying 
spin effects, an n-type wide modulation-doped quantum 
well. Due to the attraction to the ionized donors in the 
barriers, we obtain two weakly interacting electron gases 
mainly localized to the interface regions. The Rashba 
effect has often been studied experimentally and theo- 
retically in structures where the electric field is uniform 
or close to uniform. In wide modulation-doped quan- 
tum wells the electric field is strongly nonuniform [see 
Fig. [Ha)] • I n each interface region there is a strong elec- 
tric field which, according to Poisson's law, is propor- 
tional to the charge transferred from the donors in the 
barrier to the quantum well. This field is capable of pro- 
ducing a substantial Rashba effect. 

The modulation-doped quantum well is a very versatile 
system. The degree of interaction between the electron 
gases can be controlled by the well width and the carrier 
concentration. The asymmetry can be regulated by an 
applied bias or by choosing unequal spacer layer widths. 

In this paper we consider an 80 nm wide n-type 
Ino.74Gao.26Sb quantum well with Ino.7Alo.3Sb barriers. 
The growth direction is [001]. The input parameters are 
given in Table I. Among the common III-V semiconduc- 
tors InSb has the strongest spin-orbit coupling. By mix- 
ing in Ga and Al, respectively, we obtain a quantum well 
with an almost equally strong spin-orbit coupling. This 
InGaSb/InAlSb system has recently been studied exper- 
imentally by Akabori et al^ Similar effects should occur 
also for smaller well widths but the effects we want to 
display become quite clear for a well width of 80 nm. 
We here make a much more thorough analysis than in a 
recent preliminary paper M- 

We use the well established approach with an 8 x 8 
matrix approach with a minor modification (see Table 
Caption I) of the matrix given as Table C.8 in Ref. We 
apply the Burt-Foreman boundary conditions ) 12 ' 13 ! 14 and 
to avoid spurious solutions, we use a quadrature method 
in which unphysically large k values do not enter ^ The 
subband problem is solved self-consistcntly in the Hartree 
approximation. 

For reference we first consider a symmetric quantum 
well. If we, for the moment, ignore the Dresselhaus effect, 
we have a ground state with a symmetric wave function 
and a small energy separation to an excited state with an 
antisymmetric wave function. If the modulation-doped 
quantum well is wide enough that the electron gases can 
be considered as noninteracting, we seem to have a para- 
dox. Each electron gas is in a strongly asymmetric po- 
tential and a strong Rashba effect can be expected. On 



TABLE I: Input parameters used in the present work for an 
Ino.74Gao.26Sb quantum well surrounded by Ino.7Alo.3Sb bar- 
riers. We give the "true" Luttinger parameters although mod- 
ified Luttinger parameters are included in the 8x8 matrix. 1 
The conduction-band offset AE C and valence-band offset AE V 
between the two materials are also given. Mass parameters 
are given in units of the free-electron mass. Energies are 
given in electron volts. The coefficients C and Bg v describing 
the Dresselhaus effect are given in eV A and eV A 2 , respec- 
tively. We use the reasonable approximations Bj v ~ B% v and 
Bg v ps that are justified by Table 6.3 in Ref|l|. Finally we 
give the static dielectric constant e. 



Parameter InGaSb 




InAlSb 


m e 


0.01595 




0.04221 


7i 


23.87 




11.652 


72 


10.01 




4.110 


73 


11.34 




5.388 


E p 


24.27 




21.92 


E B 


0.3068 




0.790 


A 


0.777 




0.717 


C 


-0.006715 




-0.006524 


Bl 


20.596 




7.21 


AE C 




0.3344 




AE V 




0.1487 




e 


17.17 




16.00 



the other hand, if the quantum well is considered as a 
whole, the potential is symmetric and a twofold spin de- 
generacy should result [see Fig. UKa)]. This was sorted 
out previously 39 for a p-type quantum well, where the 
Rashba effect can be made several orders of magnitude 
larger than for electrons (cf. Ref. I33T ). However, the 
qualitative features are the same for the n-type quantum 
well. We recapitulate here the essentials and refer to Ref. 
l39i for details. 

First it should be noted that the signs of the electric 
fields at the interfaces are opposite to each other. If we 
label the upper spin subband at the left interface by spin 
up, the corresponding spin subband at the right inter- 
face should be labeled spin down. Looking now at the 
whole quantum well with two electron gases, the lower 
spin subband has equal amounts of spin up and spin 
down, and the twofold spin degeneracy expected for a 
symmetric potential prevails [lower subband structure in 
Fig. da)]. This implies that what looks like a spin-split 
subband in the single interface case should actually be 
considered here as two separate subbands. At k — we 
have a four-fold degeneracy because we have two spin 
directions and two electron gases. If we now consider 
a narrower but still symmetric quantum well such that 
the electron gases start to interact, the main and some- 
what unexpected effect is that the degeneracy at k = is 
partially lifted, and for all fc-values we have two closely 
spaced subbands, each with a twofold spin degeneracy 
[Fig. lib)]. 

This effect can possibly be related to the splitting re- 
cently obtained by Bernardes et who derived some 
off-diagonal matrix elements between wave functions of 
even and odd parity. This gives rise to a splitting qual- 
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0.10 



(b) 



(c) 





FIG. 1: (Color online) Schematic subband dispersions and 
conduction-band profiles for n-type MD quantum wells. The 
numbers in parentheses denote the degeneracy of the sub- 
bands for k 7^ and k = 0. In these figures the Dresselhaus 
splitting has been ignored, (a) Here the well is wide enough 
that the two electron gases at the interfaces can be consid- 
ered as noninteracting. To the left and the right we display 
the subband dispersions when the interface regions are con- 
sidered separately. The middle figure shows the case when the 
quantum well is considered as a whole, (b) Symmetric MD 
quantum well of intermediate width. The previous fourfold 
degeneracy at k = partially lifted due to the interaction 
between the electron gases, (c) Asymmetric MD quantum 
well. Due to the overall asymmetry the degeneracy is lifted 
for k / 0. 



itatively similar to the splitting described above when 
the electron (or hole) gases start to interact. However, 
the mechanism for the splitting derived in Ref. [4^ is the 




k y (1/nm) 

FIG. 2: (a) Energy subband splitting for an 80 nm n-type 
modulation-doped symmetric InGaSb quantum well. The 
wave vector is in the [10] direction in the two-dimensional 
Brillouin zone, [(b) and (c)] k dependence of the x and y 
components, respectively, of the spin vector. Here both the 
spacer layer widths are 45 nm and the electron density is 
6.8 x 10" cm' 2 . 



spin-orbit interaction. Since well widths and other input 
parameters are not given in this paper, we cannot inves- 
tigate if this effect is implicitly included in our multiband 
approach. 

The two-fold spin degeneracy becomes lifted for finite 
k if the potential of the quantum well is made asymmet- 
ric [Fig. [He)]. This is the case treated in this paper. 
The behavior described above is confirmed by numerical 
calculations.— 

We next include the Dresselhaus terms in the matrix 
but keep the quantum well potential symmetric. The 
Dresselhaus effect becomes rather different in a wide 
modulation-doped quantum well compared to a square 
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TABLE II: Spin splittings for different electric fields in various approximations. e a „ e is the voltage across the quantum well 
divided by the well width and AE ave is the energy spin splitting obtained by inserting this electric field into Eq. (6). (AE)1 xcl 
are the results when the expectation value of the electrostatic field in the layers excluding interface contributions (see text) 
averaged over filled states for spin subband i is inserted into Eq. (6). {AE)i are the corresponding results with inclusion of 
interface contributions from the valence band offsets (Eq.(7)). AE„ U m is our numerical result with the spin splittings evaluated 
at the Fermi wave vector fci? ~ 0.f3 nm _1 . The last row corresponds to the situation in Fig. 3. 



e ave AEave {AE)H cl {AE)^ cl <A£) U {AE) 11 AE nur 
(kV/cm) (meV) (meV) (meV) (meV) (meV) (meV) 



0.297 


0.012 


0.862 


-0.292 


1.384 


-0.495 


0.720 


1.242 


0.052 


1.112 


0.968 


1.768 


1.491 


2.158 


3.012 


0.135 


1.357 


1.328 


2.145 


2.040 


2.498 


4.137 


0.190 


1.485 


1.459 


2.336 


2.231 


2.646 



weliji In Fig. [U^a) we display the energy spin splitting as 
a function of wave vector. It is seen that it first rapidly 
increases but then decreases and for larger wave vectors 
(approximately above the Fermi wave vector fcj?), it stays 
rather constant at a low value. Our results imply that 
the change in spin splitting between a symmetric and an 
asymmetric wide modulation-doped quantum well nor- 
mally is dominated by the Rashba effect. In Figs. [2]^b) 
and^c) we show the k dependence of the x and y com- 
ponents of the expectation value of the spin vector X The 
absolute value of the x component decreases rapidly and 
becomes small above kp. This is in contrast to a square 
well where it stays constant^ The y component has the 
reverse behavior: small for small k values and increases 
rapidly near kp. Thus the spin direction changes from 
the x direction to the y direction as k increases along the 
[10] direction. For sufficiently large k we find a localiza- 
tion of the wave functions to one of the interface regions, 
similarly to what is shown below for the Rashba effect. 

In a quantum well with small asymmetry, the Rashba 
splitting can be comparable to the splitting between the 
symmetric and the antisymmetric states. In addition to 
a gradual transfer of wave-function amplitude to one of 
the interface regions, there will be interesting anticrossing 
phenomena, especially when k is in the [11] directionj 42 i 43 
In Table II we show the spin splittings for some values 
of the bias in various cases. Especially for small bias the 
expectation values can be quite different and even have 
opposite signs for the two spin subbands. This presents 
a fundamental problem in applying Eqs. (5) and (6) for 
more complex situations. Many downfolding approaches 
lead to expectation values of the electric field with respect 
to a subband rather than a spin subband. If one wants to 
treat this effect perturbatively, it would be appropriate to 
apply degenerate perturbation theory. However, since we 
find that a small perturbation can give rise to a substan- 
tial effect on the wave functions, the use of perturbation 
theory appears dubious in the present case. The anti- 
crossing phenomena including the contrasting results for 
k in the [11] direction have been examined more closely 
elsewhere! 42 ! 43 

In this paper we focus on the case with a fairly small 
but sufficiently large asymmetry that each wave function 
becomes almost completely localized to one of the inter- 
face regions. In the present case a bias over the quantum 
well [henceforth denoted quantum well bias (QWB)] of 



33 mV is sufficient to reach this situation. It corresponds 
to a rather small average electric field of 4.1 kV/cm (last 
row in Table II). It gives an energy separation of 9.2 
meV at k — between the lowest and next lowest sub- 
band pairs. It is shown in Fig. [3fb) that it gives a spin 
splitting that is an order of magnitude larger than for the 
same uniform electric field. By comparing columns 2 and 
7 in Table II, we can observe for different biases how the 
spin splitting is enhanced in modulation-doped quantum 
wells compared to undoped quantum wells with the same 
QWB. For the largest (last row) and the smallest biases 
(first row) , we have an enhancement by a factor of 14 and 
60, respectively. 

We thus have a modified and very efficient mechanism 
to apply a moderate QWB, and take advantage of the 
much stronger built-in electric field to obtain a substan- 
tial Rashba splitting. A qualitative explanation of the 
enhancement can be seen in Fig. 02a) where the ground- 
state wave function has become localized to one of the 
interface regions. There the electric field becomes quite 
strong and it is this local field that determines the size 
of the spin splitting. 

We note in Table II than we can reach a spin splitting 
of 2.6 meV at a wave vector of 0.13 nm" 1 . According to 
Eq. (5) this corresponds to value as large as 200 meV A 
for the Rashba parameter a. This compares well in com- 
parison to experimentally determined a valueSi 20 ! 28 ' 29 In 
addition, it is essential that in our case this a value is 
reached with a moderate bias. 

We will now investigate if we can reproduce this strong 
enhancement using the common Rashba model together 
with Eq. (6). As a first step we insert into Eq. (6) the 
expectation value of the electric field in the well and bar- 
riers ignoring any interface contributions (columns 3 and 
4). It can be expected to be enhanced by the localiza- 
tion of the wave function. In the last row of Table II the 
spin subbands have almost the same expectation value. 
However, this procedure gives a clearly smaller Rashba 
splitting than obtained in our numerical calculations. We 
have found that inclusion of BIA has a small effect on the 
results in Table II. 

As a next step it is natural to examine if the discrep- 
ancy can be explained by interface contributions. Af- 
ter solving the multiband problem in our approach, it is 
straightforward to evaluate the expectation value of the 
electric field in the layers and that of the steps at the 
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interfaces separately. If the Rashba model Eqs. (5) and 
(6) is to be used, the most reasonable approach seems to 
be to follow Ref. H3 and take the expectation value of the 
valence-band profile. The interface contributions involve 
the derivative of the discontinuities in the valence bands 
and yield Dirac delta functions in the integrals that nev- 
ertheless can be evaluated. The result becomes 

< e >interfaces= AE v (\iP(-a)\ 2 - \lf) (a) | 2 ) /e, (7) 

where the interfaces are taken at z = ±a. It is seen in Ta- 
ble II that the expectation values are increased by about 
50% when the contribution from the interface steps is 
added. This is compatible with the analytical results by 
Yang and Chang22, who concluded that the interface con- 
tribution to the Rashba splitting was typically a factor of 
two smaller than that of the electric field in the quantum 
well. 

We have also calculated the expectation value of the 
conduction-band profile to obtain the "average electric 
field" The contribution from the conduction-band off- 
sets is of opposite sign compared to that from the electric 
field in the quantum well and results in a net contribu- 
tion of 10%-15% of our numerical value. In Ref. [26| a 
contribution of only 3% was obtained. 

To examine the interface effect we have also replaced 
the barrier material by Ini-^Al^Sb of different compo- 
sitions. This changes the conduction- and valence-band 
offsets, the effective electron mass, and the nonparabol- 
icity. If the interface contribution were very sensitive to 
these parameters, this should give a clear effect. How- 
ever, in our calculations we only find a minor change. 
For example, if we change the barrier material from 
Ino.7Alo.3Sb to Ino.5Alo.5Sb (ignoring introduced strain), 
we obtain a change of the spin splitting by about 1%. 

The failure to reproduce our numerical results using 
Eqs. (5) and (6) and various expectation values raises the 
question if more elaborate downfolding approaches can 
yield better agreement or if approaches using expecta- 
tion values simply are insufficient. A more comprehensive 
comparison between various multiband and downfolding 
approaches will be published elsewhere. It is conceivable 
that two potentials having the same expectation value of 
the electric field can yield different Rashba spin splittings 
and that the spatial variation in the electric field must 
be taken into account in a multiband approach rather 
than basing the calculations on some kind of expectation 
values. 



IV. IMPLICATIONS FOR SPINTRONIC 
DEVICES 

The strong enhancement of the Rashba splitting de- 
scribed in Fig. [3] due to modulation doping can be ex- 
pected to have important implications for several spin- 
tronic devices based on the Rashba effect. For the mo- 
ment we focus on one of the best known spintronic de- 
vices, the spin transistor proposed by Datta and Das,— 
including proposed modifications of it. We will return to 
the problems encountered to make it function and first 
address the question: If it can be made to function, does 
it have the potential to become competitive with state- 
of-the-art conventional transistors? Then it is not only 



> 
LU 



0.6 



0.5 



0.4 



0.3 



1 1 

. (a) 
J 


- ' 1 1 

k 

A V " : 


- 1 1 1 — 




Y \ 

l \> 

1 





-100 -50 50 100 

z (nm) 



288.8 



288.7 



E 

LU 



288.6 - 



288.5 




4 6 
k..(1/um) 



FIG. 3: (Color online) (a) Potential, squared wave function 
and charge density, and (b) subband dispersion along the [10] 
direction for the lowest subband pair in an 80 nm InGaSb 
quantum well. The quantum well bias (potential difference 
between the interfaces) is 33 mV. Dashed lines: uniform elec- 
tric field; solid lines: modulation-doped quantum well with 
an electron density of 6.8 x 10 11 cm -2 . For the latter case the 
dotted line shows the charge distribution. 



essential that one can achieve a large wave vector split- 
ting Ak of a spin-split subband but also that it can be 
done with a small bias. As a benchmark for the per- 
formance, we choose the switch energy for Si MOSFETs 
where 3 aJ has been projected^ 

We have previously^ approximated the switch energy 
for n-type and p-type spin transistors by CV 2 , where C 
is the capacitance of a QW structure surrounded by two 
gates and V is the applied bias between them. (We have 
included here turning on and off of the device, which 
cancels a factor 1/2). We then concluded that n-type 
spin transistors with the original design would have prob- 
lems in becoming competitive with conventional transis- 
tors unless fundamentally new ideas were presented. 

A similar conclusion was independently drawn by 
Bandyopadhyay and Cahay4^ They assumed that a spin 
transistor must be based on a one-dimensional channel 
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and that only the lowest one-dimensional subband should 
be filled. However, this resulted in an anomalously small 
carrier density, 3 x 10 10 cm~ 2 or 3 x 10 5 cm" 1 . This made 
their comparison very unfavorable for the spin transistor. 

A more recent comparison with conventional transis- 
tors has been made by Hall and Flatte 47 for a modi- 
fied spin transistor. It is not based on the Rashba ef- 
fect but rather on gate-induced spin relaxation. A cru- 
cial factor in their approach seems to be efficiency of 
the gate-induced spin relaxation compared to other spin- 
relaxation mechanisms. Their comparison was quite fa- 
vorable for the spin transistor. They estimated a switch 
energy of 0.5 aJ which is similar to what we find be- 
low for our modified spin transistor. The performance 
of this transistor has been subject to some controversy 
concerning the need for very efficient spin injection. 4 ® 
Since the Hall-Flatte spin transistor is based on another 
mechanism than the spin precession considered here, it 
is beyond the scope of the present paper to enter this 
debate. 

Utilizing the built-in electric field in the modulation- 
doped quantum well, one can achieve a given Afc with a 
QWB that is an order of magnitude smaller than with 
a uniform electric field. If we only consider the lowest 
spin subband pair and follow the approach of Ref. |45| . we 
obtain a switch energy of 0.4 aJ in the modulation-doped 
case and 35 aJ in a spin transistor with the same length 
and uniform electric field. The former figure compares 
very well with present state-of-the-art transistors. Thus 
the utilization of the mechanism proposed by us could 
make a substantial difference for the competitiveness of 
spin transistors. 

We have calculated the additional contribution to the 
switch energy from the redistribution of carriers in the 
QW, taking the fey -dependent wave functions into ac- 
count but found that it only increases by about 20 %. 

A complication with our design is that the second sub- 
band pair with the opposite sign of Afc and spin preces- 
sion direction is also filled. This does not prevent the 
possibility that the spins at the two interfaces can have 
made a precession by the angle it but in opposite direc- 
tions on the arrival to the drain where the transmission 
becomes low. 

It has been demonstrated that one can contact the 
electron gases in a double quantum well structure 
separately.— It seems feasible that also the interface 
regions of a wide modulation-doped quantum well can 
be contacted separately which opens up interesting pos- 
sibilities occurring from the controllable properties of 
modulation-doped quantum wells. 

One can envision practical problems to create a per- 
fectly symmetric quantum well structure corresponding 
to the on state of a spin transistor. One possibility is a 
double-gate structure in which the total carrier concen- 
tration and the asymmetry can be controlled separately. 
In Ref. |50| the back gate voltage was of the order 100 V, 
which is not very practical for devices. An alternative 
design 45 is to have a heavily doped semiconductor layer 
just below the quantum well structure. In this way a 
larger fraction of the applied voltage falls over the quan- 
tum well. 

We now turn to the problem of making a spin transistor 



function, possibly with some modification of the original 
idea.— A fundamental problem is that the Rashba effect 
can be described in terms of an effective magnetic field 
that is perpendicular to both the electric field and the di- 
rection of motion for the carrier #i Even spin-independent 
scattering leads to a change of the direction of the ve- 
locity and thus the axis of the spin precession. It can 
also be difficult to inject all the carriers in the same di- 
rection. In the case we consider in Fig. [3l we obtain a 
precession length L = ir/Ak ~ 1/im. Ballistic trans- 
port over such a distance requires rather low tempera- 
tures. An idea with the purpose of balancing the Rashba 
and Dresselhaus effect 5 -^ by setting a = (3 in the lin- 
earized model makes diffusive transport possible but at 
the price of a substantial transmission in the off state. 
One-dimensional channels have been proposed in which 
the carriers are more or less forced to move in the same 
direction. As mentioned above the small energy sepa- 
ration between the one-dimensional subbands leads to 
multi-mode transport for realistic carrier densities. This 
is not a prohibitive problem as has been demonstrated 
by Lusakowski et al£& 

A fundamental problem is that in the approach with 
the Rashba term, which is a reasonable approximation 
for electrons in an undoped quantum well, the Rashba 
splitting Afc is proportional to the Rashba coefficient a 
but the spin decoherence rate becomes proportional to 
a (Ref. 3). A large a is beneficial for a rapid spin 
precession and corresponds to a short gate length in a 
spin transistor but this advantage is thus offset by the 
shorter spin decoherence time. 

An alternate approach has been presented by Bandy- 
opadhyay and Cahay. 53 Instead of relying on the Rashba 
effect, they propose using the Dresselhaus effect in a 
structure with a split gate and a parabolic potential. The 
main reason was to avoid an in-plane magnetic field in 
the semiconductor from the magnetized source and drain. 
However, the main requirement for spin precession is that 
their magnetization is perpendicular to the effective mag- 
netic field in the channel. Thus it can be either along 
the channel (as drawn in Ref. IBH ) or perpendicular to 
the layers (as drawn in Ref. I4lh . The mechanism in the 
transistor based on the Dresselhaus effect is changing the 
bias of the split gate and then it is assumed that the 
curvature of the parabolic potential changes. Numerical 
calculations 5 - 4 - have indicated that, however, the effect of 
changing the bias is mainly that the potential becomes 
flatter in the middle when the channel starts to fill while 
the curvature of the side walls does not change much. 

An alternative that has not been given much atten- 
tion so far is a p-type spin transistor. With a suit- 
able design we have shown that one can obtain a large 
Afc with an electric field as small as 2 kV/cm (Ref. 
[33l ). The corresponding precession length is only 40 
nm and the possibility of having ballistic transport over 
such a short distance clearly seems feasible. The strong 
spin-orbit interaction including its dependence on a gate 
voltage has recently been demonstrated for a p-type 
GaAs/AlGaAs heterostructure by Grbic et a/.— Accord- 
ing to our calculations^ 3 - even stronger spin splitting can 
be achieved for higher hole densities. Furthermore, the 
strong anisotropy of hole subbands can possibly be uti- 
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lized to get a preferred direction of motion without lat- 
eral confinement. For holes there is no simple relation 
between spin precession and spin decoherence rates. Es- 
timates based on experimental determinations indicate 
that the spin decoherence time can be much longer than 
the transit time. Because of the strongly nonparabolic 
hole subbands and their mixed heavy-hole and light-hole 
character, rather cumbersome numerical calculations ap- 
pear necessary for a more accurate prediction of the 
transport properties. For small k analytical expressions 
proportional to k 3 for the Rashba splitting in heavy-hole 
subbands have been derived* 1 - However, it has been found 
that the largest spin splittings occur beyond the range of 
validity of this models 

A relevant question is if one can combine the superef- 
ficient Rashba effect for holes with the enhancement in 
modulation-doped quantum wells presented here. How- 
ever, we have shown that for p-type spin transistors the 
largest Rashba splitting is obtained for quite small elec- 
tric fields (~ 5 - 10 kV/cm) while the effect of modula- 
tion doping is to apply a small bias to utilize the built-in 
electric field of the order 50 — 100 kV/cm for which the 
Rashba effect for holes is reduced. 

A well-known problem is that the conductivity mis- 
match between metal and semiconductor severely limits 
the spin injection efficiency— A proposal by Rashba is 
having tunnel barriers between the metal contacts and 
the semiconductor^ A fundamental problem recently 
pointed out by Fert et al— is that this decreases the 
transmission coefficient and increases the dwell time such 
that it can become long compared to the spin dephasing 
time in semiconductor-based spin transistors. They in- 
stead proposed using carbon nanotubes. It is beyond the 
scope of the present paper to evaluate the competitive- 
ness of semiconductors vs. carbon nanotubes for spin- 
tronic devices. However, we would like to point out that 
this problem occurs for injection from a spin-polarized 
contact but other solutions in the form of spin filters 
have been proposed. A particularly interesting idea is 
to put a magnetic layer on top of a layered semiconduc- 
tor structure such that the in-plane fringe fields act as 
a spin filter^ The appealing aspect of this solution is 
that current flows in the channel below the metal with- 
out passing any interfaces where the spin polarization can 
be reduced. 



V. DISCUSSION AND CONCLUSIONS 

We have implicitly assumed coherence of the wave 
function across the 80 nm QW with a high and broad 
barrier in the middle. Whether this coherence actually 
prevails should depend on the sample quality. This sys- 
tem with our predicted effects seems ideal for further 
studies of this fundamental problem. 

The self-consistent calculations have so far been per- 
formed in the Hartree approximation. For studies of spin 
properties it is conceivable that exchange and correlation 
can give significant effects, especially in anticrossing sit- 
uations. This is planned to be examined in future publi- 
cations. 

In our multiband approach the well established Burt- 



Foreman boundary condition o 12 i 13 ' are behind any in- 
terface contribution. The exact relation behind this ap- 
proach and what is obtained by folding down the adja- 
cent bands to the conduction band as in Refs. [26| and 
l27l . is not trivial and remains to be analyzed. Our ap- 
proach is in our opinion more sound than one-band ap- 
proaches that are based on approximations whose accu- 
racy is difficult to determine. Furthermore, the analytical 
expression s] 26 ! 27 are based on consideration of the inter- 
faces, where the EFA has its main weakness and where 
the actual gradual but rapid potential variation near an 
interface is replaced by a sharp step. Especially when the 
interfaces give substantial contributions, it is likely that 
different operator orderings can influence the results con- 
siderably. Calculations with downfolding procedures for 
wide modulation-doped quantum wells and comparison 
with the present results would be valuable in evaluating 
how close the results of these approaches are in a non- 
trivial case such as this. 

It has recently been predicted that the spin Hall ef- 
fect can be strongly enhanced at a subband anticrossing 
in a bilayer system.— There the potential was not spec- 
ified but the Rashba coefficients were allowed to differ 
in the two layers. For further investigations of this ef- 
fect modulation-doped quantum wells seem useful due to 
the possibility of controlling the degree of interaction be- 
tween the two electron gases and each of the interface 
fields. 

In conclusion we have analyzed the foundations of 
the envelope function approximation and concluded that, 
while the Dresselhaus effect should be included as k- 
dependent terms in a matrix, the proper inclusion of the 
Rashba effect is adding the macroscopic potential along 
the diagonal in a multiband approach. This has given 
good agreement with experiment for two-dimensional 
hole gases . 21 ' 23 i 24 The commonly used Rashba term (5) 
is a hybrid including both potential and k. The proper 
derivation of such a term within the framework of the 
EFA with proper boundary condition o 12 ^ 13 ' 14 deserves to 
be examined more closely. 

For symmetric wells with Dresselhaus effect only, we 
find interesting effects in a modulation-doped quantum 
well that are qualitatively different from those in a square 
well. 

We have found that, with a non- uniform electric field, 
insertion of some kind of expectation value or other av- 
erage into Eq. (6) underestimates the Rashba splitting. 
Furthermore, this expectation value is not always well- 
defined for a subband because it can differ substantially 
between its spin-split components. We have found that 
the contribution from the interfaces is about half of that 
from the electric field in the layers for the potential we 
have considered. 

We have demonstrated a very efficient switching mech- 
anism of the Rashba splitting in wide modulation-doped 
quantum wells. One can use a bias corresponding to 
a moderate average electric field and still get a Rashba 
splitting typically enhanced by an order of magnitude due 
to the built-in local electric field in the interface region. 
The switching mechanism is based on localization of each 
wave function to one interface region with a barely suffi- 
cient bias. A switching mechanism based on anticrossing 
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in slightly asymmetric quantum well o i 3 is not included 
here but will be examined further elsewhere. 

The enhancement of the efficiency of the Rashba 
effect should be valuable for different spintronic devices. 
We here have focused on spin transistors of the type 
proposed by Datta and Das.- With our modifica- 
tion we find that it can get a potential to outperform 
conventional transistors. We have also discussed some re- 



maining obstacles to make such spin transistors function. 
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